HIV-1 capsid shape, orientation, and entropic elasticity regulate translocation into the nuclear pore complex.

Nuclear import and uncoating of the viral capsid are critical steps in the HIV-1 life cycle that serve to transport and release genomic material into the nucleus. Viral core import involves translocating the HIV-1 capsid at the nuclear pore complex (NPC). Notably, the central channel of the NPC appears to often accommodate and allow passage of intact HIV-1 capsid, though mechanistic details of the process remain to be fully understood. Here, we investigate the molecular interactions that operate in concert between the HIV-1 capsid and the NPC that regulate capsid translocation through the central channel. To this end, we develop a "bottom-up" coarse-grained (CG) model of the human NPC from recently released cryo-electron tomography structure and then construct composite membrane-embedded CG NPC models. We find that successful translocation from the cytoplasmic side to the NPC central channel is contingent on the compatibility of the capsid morphology and channel dimension and the proper orientation of the capsid approach to the channel from the cytoplasmic side. The translocation dynamics is driven by maximizing the contacts between phenylalanine-glycine nucleoporins at the central channel and the capsid. For the docked intact capsids, structural analysis reveals correlated striated patterns of lattice disorder likely related to the intrinsic capsid elasticity. Uncondensed genomic material inside the docked capsid augments the overall lattice disorder of the capsid. Our results suggest that the intrinsic "elasticity" can also aid the capsid to adapt to the stress and remain structurally intact during translocation.

Dynamics of upstream ESCRT organization at the HIV-1 budding site.

In the late stages of the HIV-1 life cycle, membrane localization and self-assembly of Gag polyproteins induce membrane deformation and budding. Release of the virion requires direct interaction between immature Gag lattice and upstream ESCRT machinery at the viral budding site, followed by assembly of downstream ESCRT-III factors, culminating in membrane scission. However, molecular details of upstream ESCRT assembly dynamics at the viral budding site remain unclear. In this work, using coarse-grained (CG) molecular dynamics (MD) simulations, we investigated the interactions between Gag, ESCRT-I, ESCRT-II, and membrane to delineate the dynamical mechanisms by which upstream ESCRTs assemble templated by late-stage immature Gag lattice. We first systematically derived "bottom-up" CG molecular models and interactions of upstream ESCRT proteins from experimental structural data and extensive all-atom MD simulations. Using these molecular models, we performed CG MD simulations of ESCRT-I oligomerization and ESCRT-I/II supercomplex formation at the neck of the budding virion. Our simulations demonstrate that ESCRT-I can effectively oligomerize to higher-order complexes templated by the immature Gag lattice both in the absence of ESCRT-II and when multiple copies of ESCRT-II are localized at the bud neck. The ESCRT-I/II supercomplexes formed in our simulations exhibit predominantly columnar structures, which has important implications for the nucleation pathway of downstream ESCRT-III polymers. Importantly, ESCRT-I/II supercomplexes bound to Gag initiate membrane neck constriction by pulling the inner edge of the bud neck closer to the ESCRT-I headpiece ring. Our findings serve to elucidate a network of interactions between upstream ESCRT machinery, immature Gag lattice, and membrane neck that regulate protein assembly dynamics at the HIV-1 budding site.

Prion-like low complexity regions enable avid virus-host interactions during HIV-1 infection.

Cellular proteins CPSF6, NUP153 and SEC24C play crucial roles in HIV-1 infection. While weak interactions of short phenylalanine-glycine (FG) containing peptides with isolated capsid hexamers have been characterized, how these cellular factors functionally engage with biologically relevant mature HIV-1 capsid lattices is unknown. Here we show that prion-like low complexity regions (LCRs) enable avid CPSF6, NUP153 and SEC24C binding to capsid lattices. Structural studies revealed that multivalent CPSF6 assembly is mediated by LCR-LCR interactions, which are templated by binding of CPSF6 FG peptides to a subset of hydrophobic capsid pockets positioned along adjoining hexamers. In infected cells, avid CPSF6 LCR-mediated binding to HIV-1 cores is essential for functional virus-host interactions. The investigational drug lenacapavir accesses unoccupied hydrophobic pockets in the complex to potently impair HIV-1 inside the nucleus without displacing the tightly bound cellular cofactor from virus cores. These results establish previously undescribed mechanisms of virus-host interactions and antiviral action.

A helical assembly of human ESCRT-I scaffolds reverse-topology membrane scission.

The ESCRT complexes drive membrane scission in HIV-1 release, autophagosome closure, multivesicular body biogenesis, cytokinesis, and other cell processes. ESCRT-I is the most upstream complex and bridges the system to HIV-1 Gag in virus release. The crystal structure of the headpiece of human ESCRT-I comprising TSG101-VPS28-VPS37B-MVB12A was determined, revealing an ESCRT-I helical assembly with a 12-molecule repeat. Electron microscopy confirmed that ESCRT-I subcomplexes form helical filaments in solution. Mutation of VPS28 helical interface residues blocks filament formation in vitro and autophagosome closure and HIV-1 release in human cells. Coarse-grained (CG) simulations of ESCRT assembly at HIV-1 budding sites suggest that formation of a 12-membered ring of ESCRT-I molecules is a geometry-dependent checkpoint during late stages of Gag assembly and HIV-1 budding and templates ESCRT-III assembly for membrane scission. These data show that ESCRT-I is not merely a bridging adaptor; it has an essential scaffolding and mechanical role in its own right.

The end of ice I.

The appearance of ice I in the smallest possible clusters and the nature of its phase coexistence with liquid water could not thus far be unraveled. The experimental and theoretical infrared spectroscopic and free-energy results of this work show the emergence of the characteristic hydrogen-bonding pattern of ice I in clusters containing only around 90 water molecules. The onset of crystallization is accompanied by an increase of surface oscillator intensity with decreasing surface-to-volume ratio, a spectral indicator of nanoscale crystallinity of water. In the size range from 90 to 150 water molecules, we observe mixtures of largely crystalline and purely amorphous clusters. Our analysis suggests that the liquid-ice I transition in clusters loses its sharp 1st-order character at the end of the crystalline-size regime and occurs over a range of temperatures through heterophasic oscillations in time, a process without analog in bulk water.

Hydrogen-Bonding and Hydrophobic Groups Contribute Equally to the Binding of Hyperactive Antifreeze and Ice-Nucleating Proteins to Ice.

Hyperactive insect antifreeze proteins and bacterial ice-nucleating proteins are arguably the most potent ice-binding molecules in nature. These highly effective proteins bind ice through amphiphilic ice-binding sites based on arrays of threonine residues. It remains poorly understood how hydrophilic and hydrophobic groups of the binding site contribute to the ice affinity of proteins. Here, we use molecular simulations to demonstrate that the hydrogen-bonding and hydrophobic groups at the ice-binding site of the antifreeze protein TmAFP of Tenebrio molitor and extended ice-nucleating protein surfaces contribute distinctively yet almost equally in magnitude to their binding free energy to ice. The methyl groups rigidize the ice-binding site, slow the water dynamics at the ice-binding surface, and stabilize the clathrate-like water in the anchored clathrate motif that binds these proteins to ice. We find that hydrophobic dehydration of the methyl group does not contribute to the binding free energy of the protein to ice. The role of the hydroxyl groups is to anchor the clathrate-like water through direct hydrogen-bonding, positioning and slowing the dynamics of water at the trough of the binding site. We uncover a correlation between slower dynamics of water at the binding site for the protein in solution and stronger free energy of binding of the protein to ice. The synergy between hydrophobic and hydrophilic groups unveiled by this study provides guidance for the design of synthetic ice-binding molecules with enhanced ice nucleation and antifreeze activity.

How Size and Aggregation of Ice-Binding Proteins Control Their Ice Nucleation Efficiency.

Organisms that thrive at cold temperatures produce ice-binding proteins to manage the nucleation and growth of ice. Bacterial ice-nucleating proteins (INP) are typically large and form aggregates in the cell membrane, while insect hyperactive antifreeze proteins (AFP) are soluble and generally small. Experiments indicate that larger ice-binding proteins and their aggregates nucleate ice at warmer temperatures. Nevertheless, a quantitative understanding of how size and aggregation of ice-binding proteins determine the temperature Thet at which proteins nucleate ice is still lacking. Here, we address this question using molecular simulations and nucleation theory. The simulations indicate that the 2.5 nm long antifreeze protein TmAFP nucleates ice at 2 ± 1 °C above the homogeneous nucleation temperature, in good agreement with recent experiments. We predict that the addition of ice-binding loops to TmAFP increases Thet, but not enough to compete in efficiency with the bacterial INP. We implement an accurate procedure to determine Thet of surfaces of finite size using classical nucleation theory, and, after validating the theory against Thet of the proteins in molecular simulations, we use it to predict Thet of the INP of Ps. syringae as a function of the length and number of proteins in the aggregates. We conclude that assemblies with at most 34 INP already reach the Thet = -2 °C characteristic of this bacterium. Interestingly, we find that Thet is a strongly varying nonmonotonic function of the distance between proteins in the aggregates. This indicates that, to achieve maximum freezing efficiency, bacteria must exert exquisite, subangstrom control of the distance between INP in their membrane.

Preordering of water is not needed for ice recognition by hyperactive antifreeze proteins.

Antifreeze proteins (AFPs) inhibit ice growth in organisms living in cold environments. Hyperactive insect AFPs are particularly effective, binding ice through "anchored clathrate" motifs. It has been hypothesized that the binding of hyperactive AFPs to ice is facilitated by preordering of water at the ice-binding site (IBS) of the protein in solution. The antifreeze protein TmAFP displays the best matching of its binding site to ice, making it the optimal candidate to develop ice-like order in solution. Here we use multiresolution simulations to unravel the mechanism by which TmAFP recognizes and binds ice. We find that water at the IBS of the antifreeze protein in solution does not acquire ice-like or anchored clathrate-like order. Ice recognition occurs by slow diffusion of the protein to achieve the proper orientation with respect to the ice surface, followed by fast collective organization of the hydration water at the IBS to form an anchored clathrate motif that latches the protein to the ice surface. The simulations suggest that anchored clathrate order could develop on the large ice-binding surfaces of aggregates of ice-nucleating proteins (INP). We compute the infrared and Raman spectra of water in the anchored clathrate motif. The signatures of the OH stretch of water in the anchored clathrate motif can be distinguished from those of bulk liquid in the Raman spectra, but not in the infrared spectra. We thus suggest that Raman spectroscopy may be used to probe the anchored clathrate order at the ice-binding surface of INP aggregates.

Ice-Nucleating and Antifreeze Proteins Recognize Ice through a Diversity of Anchored Clathrate and Ice-like Motifs.

Cold-adapted organisms produce antifreeze and ice-nucleating proteins to prevent and promote ice formation. The crystal structure of hyperactive bacterial antifreeze protein (AFP) MpAFP suggests that this protein binds ice through an anchored clathrate motif. It is not known whether other hyperactive AFPs and ice-nucleating proteins (INPs) use the same motif to recognize or nucleate ice. Here we use molecular simulations to elucidate the ice-binding motifs of hyperactive insect AFPs and a model INP of Pseudomonas syringae. We find that insect AFPs recognize ice through anchored clathrate motifs distinct from that of MpAFP. By performing simulations of ice nucleation by PsINP, we identify two distinct ice-binding sites on opposite sides of the ß-helix. The ice-nucleating sequences identified in the simulations agree with those previously proposed for the closely related INP of Pseudomonas borealis based on the structure of the protein. The simulations indicate that these sites have comparable ice-nucleating efficiency, but distinct binding motifs, controlled by the amino acid sequence: one is an anchored clathrate and the other ice-like. We conclude that anchored clathrate and ice-like motifs can be equally effective for binding proteins to ice and promoting ice nucleation.

Role of stacking disorder in ice nucleation.

The freezing of water affects the processes that determine Earth's climate. Therefore, accurate weather and climate forecasts hinge on good predictions of ice nucleation rates. Such rate predictions are based on extrapolations using classical nucleation theory, which assumes that the structure of nanometre-sized ice crystallites corresponds to that of hexagonal ice, the thermodynamically stable form of bulk ice. However, simulations with various water models find that ice nucleated and grown under atmospheric temperatures is at all sizes stacking-disordered, consisting of random sequences of cubic and hexagonal ice layers. This implies that stacking-disordered ice crystallites either are more stable than hexagonal ice crystallites or form because of non-equilibrium dynamical effects. Both scenarios challenge central tenets of classical nucleation theory. Here we use rare-event sampling and free energy calculations with the mW water model to show that the entropy of mixing cubic and hexagonal layers makes stacking-disordered ice the stable phase for crystallites up to a size of at least 100,000 molecules. We find that stacking-disordered critical crystallites at 230 kelvin are about 14 kilojoules per mole of crystallite more stable than hexagonal crystallites, making their ice nucleation rates more than three orders of magnitude higher than predicted by classical nucleation theory. This effect on nucleation rates is temperature dependent, being the most pronounced at the warmest conditions, and should affect the modelling of cloud formation and ice particle numbers, which are very sensitive to the temperature dependence of ice nucleation rates. We conclude that classical nucleation theory needs to be corrected to include the dependence of the crystallization driving force on the size of the ice crystallite when interpreting and extrapolating ice nucleation rates from experimental laboratory conditions to the temperatures that occur in clouds.

Sink or Swim: Ions and Organics at the Ice-Air Interface.

The ice-air interface is an important locus of environmental chemical reactions. The structure and dynamics of the ice surface impact the uptake of trace gases and kinetics of reactions in the atmosphere and snowpack. At tropospheric temperatures, the ice surface is partially premelted. Experiments indicate that ions increase the liquidity of the ice surface but hydrophilic organics do not. However, it is not yet known the extent of the perturbation solutes induce at the ice surface and what is the role of the disordered liquid-like layer in modulating the interaction between solutes and their mobility and aggregation at the ice surface. Here we use large-scale molecular simulations to investigate the effect of ions and glyoxal, one of the most abundant oxygenated volatile organic compounds in the atmosphere, on the structure, dynamics, and solvation properties of the ice surface. We find that the premelted surface of ice has unique solvation properties, different from those of liquid water. The increase in surface liquidity resulting from the hydration of ions leads to a water-mediated attraction of ions at the ice surface. Glyoxal molecules, on the other hand, perturb only slightly the surface of ice and do not experience water-driven attraction. They nonetheless accumulate as dry agglomerates at the ice surface, driven by direct interactions between the organic molecules. The enhanced attraction and clustering of ions and organics at the ice surface may play a significant role in modulating the mechanism and rate of heterogeneous chemical reactions occurring at the surface of atmospheric ice particles.

Ice Nucleation Efficiency of Hydroxylated Organic Surfaces Is Controlled by Their Structural Fluctuations and Mismatch to Ice.

Heterogeneous nucleation of ice induced by organic materials is of fundamental importance for climate, biology, and industry. Among organic ice-nucleating surfaces, monolayers of long chain alcohols are particularly effective, while monolayers of fatty acids are significantly less so. As these monolayers expose to water hydroxyl groups with an order that resembles the one in the basal plane of ice, it was proposed that lattice matching between ice and the surface controls their ice-nucleating efficiency. Organic monolayers are soft materials and display significant fluctuations. It has been conjectured that these fluctuations assist in the nucleation of ice. Here we use molecular dynamic simulations and laboratory experiments to investigate the relationship between the structure and fluctuations of hydroxylated organic surfaces and the temperature at which they nucleate ice. We find that these surfaces order interfacial water to form domains with ice-like order that are the birthplace of ice. Both mismatch and fluctuations decrease the size of the preordered domains and monotonously decrease the ice freezing temperature. The simulations indicate that fluctuations depress the freezing efficiency of monolayers of alcohols or acids to half the value predicted from lattice mismatch alone. The model captures the experimental trend in freezing efficiencies as a function of chain length and predicts that alcohols have higher freezing efficiency than acids of the same chain length. These trends are mostly controlled by the modulation of the structural mismatch to ice. We use classical nucleation theory to show that the freezing efficiencies of the monolayers are directly related to their free energy of binding to ice. This study provides a general framework to relate the equilibrium thermodynamics of ice binding to a surface and the nonequilibrium ice freezing temperature and suggests that these could be predicted from the structure of interfacial water.

What Determines the Ice Polymorph in Clouds?

Ice crystals in the atmosphere nucleate from supercooled liquid water and grow by vapor uptake. The structure of the ice polymorph grown has strong impact on the morphology and light scattering of the ice crystals, modulates the amount of water vapor in ice clouds, and can impact the molecular uptake and reactivity of atmospheric aerosols. Experiments and molecular simulations indicate that ice nucleated and grown from deeply supercooled liquid water is metastable stacking disordered ice. The ice polymorph grown from vapor has not yet been determined. Here we use large-scale molecular simulations to determine the structure of ice that grows as a result of uptake of water vapor in the temperature range relevant to cirrus and mixed-phase clouds, elucidate the molecular mechanism of the formation of ice at the vapor interface, and compute the free energy difference between cubic and hexagonal ice interfaces with vapor. We find that vapor deposition results in growth of stacking disordered ice only under conditions of extreme supersaturation, for which a nonequilibrium liquid layer completely wets the surface of ice. Such extreme conditions have been used to produce stacking disordered frost ice in experiments and may be plausible in the summer polar mesosphere. Growth of ice from vapor at moderate supersaturations in the temperature range relevant to cirrus and mixed-phase clouds, from 200 to 260 K, produces exclusively the stable hexagonal ice polymorph. Cubic ice is disfavored with respect to hexagonal ice not only by a small penalty in the bulk free energy (3.6 ± 1.5 J mol(-1) at 260 K) but also by a large free energy penalty at the ice-vapor interface (89.7 ± 12.8 J mol(-1) at 260 K). The latter originates in higher vibrational entropy of the hexagonal-terminated ice-vapor interface. We predict that the free energy penalty against the cubic ice interface should decrease strongly with temperature, resulting in some degree of stacking disorder in ice grown from vapor in the tropical tropopause layer, and in polar stratospheric and noctilucent clouds. Our findings support and explain the evolution of the morphology of ice crystals from hexagonal to trigonal symmetry with decreasing temperature, as reported by experiments and in situ measurements in clouds. We conclude that selective growth of the elusive cubic ice polymorph by manipulation of the interfacial properties can likely be achieved at the ice-liquid interface but not at the ice-vapor interface.

Free energy contributions and structural characterization of stacking disordered ices.

Crystallization of ice from deeply supercooled water and amorphous ices - a process of fundamental importance in the atmosphere, interstellar space, and cryobiology - results in stacking disordered ices with a wide range of metastabilities with respect to hexagonal ice. The structural origin of this high variability, however, has not yet been elucidated. Here we use molecular dynamics simulations with the mW water model to characterize the structure of ice freshly grown from supercooled water at temperatures from 210 to 270 K, the thermodynamics of stacking faults, line defects, and interfaces, and to elucidate the interplay between kinetics and thermodynamics in determining the structure of ice. In agreement with experiments, the ice grown in the simulations is stacking disordered with a random distribution of cubic and hexagonal layers, and a cubicity that decreases with growth temperature. The former implies that the cubicity of ice is determined by processes at the ice/liquid interface, without memory of the structure of buried ice layers. The latter indicates that the probability of building a cubic layer at the interface decreases upon approaching the melting point of ice, which we attribute to a more efficient structural equilibration of ice at the liquid interface as the driving force for growth wanes. The free energy cost for creating a pair of cubic layers in ice is 8.0 J mol(-1) in experiments, and 9.7 ± 1.9 J mol(-1) for the mW water model. This not only validates the simulations, but also indicates that dispersion in cubicity is not sufficient to explain the large energetic variability of stacking disordered ices. We compute the free energy cost of stacking disorder, line defects, and interfaces in ice and conclude that a characterization of the density of these defects is required to predict the degree of metastability and vapor pressure of atmospheric ices.

Ice crystallization in ultrafine water-salt aerosols: nucleation, ice-solution equilibrium, and internal structure.

Atmospheric aerosols have a strong influence on Earth's climate. Elucidating the physical state and internal structure of atmospheric aqueous aerosols is essential to predict their gas and water uptake, and the locus and rate of atmospherically important heterogeneous reactions. Ultrafine aerosols with sizes between 3 and 15 nm have been detected in large numbers in the troposphere and tropopause. Nanoscopic aerosols arising from bubble bursting of natural and artificial seawater have been identified in laboratory and field experiments. The internal structure and phase state of these aerosols, however, cannot yet be determined in experiments. Here we use molecular simulations to investigate the phase behavior and internal structure of liquid, vitrified, and crystallized water-salt ultrafine aerosols with radii from 2.5 to 9.5 nm and with up to 10% moles of ions. We find that both ice crystallization and vitrification of the nanodroplets lead to demixing of pure water from the solutions. Vitrification of aqueous nanodroplets yields nanodomains of pure low-density amorphous ice in coexistence with vitrified solute rich aqueous glass. The melting temperature of ice in the aerosols decreases monotonically with an increase of solute fraction and decrease of radius. The simulations reveal that nucleation of ice occurs homogeneously at the subsurface of the water-salt nanoparticles. Subsequent ice growth yields phase-segregated, internally mixed, aerosols with two phases in equilibrium: a concentrated water-salt amorphous mixture and a spherical cap-like ice nanophase. The surface of the crystallized aerosols is heterogeneous, with ice and solution exposed to the vapor. Free energy calculations indicate that as the concentration of salt in the particles, the advance of the crystallization, or the size of the particles increase, the stability of the spherical cap structure increases with respect to the alternative structure in which a core of ice is fully surrounded by solution. We predict that micrometer-sized particles and nanoparticles have the same equilibrium internal structure. The variation of liquid-vapor surface tension with solute concentration is a key factor in determining whether a solution-embedded ice core or vapor-exposed ice cap is the equilibrium structure of the aerosols. In agreement with experiments, we predict that the structure of mixed-phase HNO3-water particles, representative of polar stratospheric clouds, consists of an ice core surrounded by freeze-concentrated solution. The results of this work are important to determine the phase state and internal structure of sea spray ultrafine aerosols and other mixed-phase particles under atmospherically relevant conditions.

Heterogeneous nucleation of ice on carbon surfaces.

Atmospheric aerosols can promote the heterogeneous nucleation of ice, impacting the radiative properties of clouds and Earth's climate. The experimental investigation of heterogeneous freezing of water droplets by carbonaceous particles reveals widespread ice freezing temperatures. It is not known which structural and chemical characteristics of soot account for the variability in ice nucleation efficiency. Here we use molecular dynamics simulations to investigate the nucleation of ice from liquid water in contact with graphitic surfaces. We find that atomically flat carbon surfaces promote heterogeneous nucleation of ice, while molecularly rough surfaces with the same hydrophobicity do not. Graphitic surfaces and other surfaces that promote ice nucleation induce layering in the interfacial water, suggesting that the order imposed by the surface on liquid water may play an important role in the heterogeneous nucleation mechanism. We investigate a large set of graphitic surfaces of various dimensions and radii of curvature and find that variations in nanostructures alone could account for the spread in the freezing temperatures of ice on soot in experiments. We conclude that a characterization of the nanostructure of soot is needed to predict its ice nucleation efficiency.